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, We propose an extension to other high T c compounds of a model introduced earlier for YBCO. 

' In the "self-doped" compounds we assume that the doping part ( namely the BiO, HgO, TIO planes 

Q\ \ in BSCCO, HBCCO, TBCCO respectively ) is metallic, which leads to a multiband model. This 

assumption is supported by band structure calculations. Taking a repulsive pairing interaction 
between these doping bands and the CuOi bands leads to opposite signs for the order parameter on 
these bands and to nodes whenever the Fermi surfaces of these bands cross. We show that in BSCCO 
the low temperature dependence of the penetration depth is reasonably accounted for. In this case 
CO ' the nodes are not located near the 45° direction, which makes the experimental determination of 

CI the node locations an important test for our model. The situation in HBCCO and TBCCO is rather 

analogous to BSCCO. We consider the indications given by NMR and find that they rather favor 
a metallic character for the doping bands. Finally we discuss the cases of NCCO and LSCO which 
Q ■ are not "self-doped" and where our model does not give nodes. 
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PACS numbers : 74.20.Fg, 74.72.Bk, 74.25.Jb 



I. INTRODUCTION 



The recent years have seen a marked progress in the debate on the mechanism of high T c superconductivity. Indeed 
' the existence of nodes and change of sign of the order parameter at the Fermi surface has found a very important 
experimental support in some compounds like Y BaiCu-sOi (YBCO) . This includes on one hand various experiments 
showing the existence of low energy excited states, such as tunneling, Raman scattering and penetration depth [Q. 
O , Moreover more spectacular experiments designed to probe the change of sign have given positive answers. These 
are the corner SQUID experiments Q and the observation of the half magnetic flux quantum in 3 grain-boundary 
y—\ \ Josephson junctions This change of sign is a clear indication that there is some important repulsive contribution 
^ • in the pairing interaction. Since the spin fluctuation mechanism Q is based uniquely on a repulsive interaction, it 
appears as a leading candidate for the pairing mechanism in these high T c superconductors. In its simplest version it 
leads to an order parameter with a d x 2_ y 2 symmetry, commonly called d-wave order parameter. 

Unfortunately this simple spin fluctuation interpretation meets some problems in YBCO. An important one is the 
^ — ■ weak sensitivity of the critical temperature of Y Ba^Cu^O'j to defects. Any kind of impurities produces in d-wave 
Q\ , superconductors an effect analogous to pair-breaking by magnetic impurities in standard s-wave superconductors, 
leading to a rapid decrease of the critical temperature with increasing scattering rate, as measured for example by the 
residual resistivity. This is in contrast with some experiments showing a decrease of T c an order of magnitude slower 
than expected M. 



Another problem is the marked anisotropy in the a-b plane shown by a number of recent experiments in YBCO. 
The existence of a sizeable Josephson current Q in a c-axis Pb — Y Ba2Cu^O-j tunnelling junction shows that pure 
q \ d-wave symmetry is not satisfied. Because of the large size of this Josephson current as compared to what would be 
expected from a superconductor with s-wave symmetry, it is difficult to believe that the departure from pure d-wave 
is only due to the small orthorhombic 1,6% distorsion of the CuC>2 planes. Various other properties have also shown 
a clear in-plane anisotropy. The resistivity has an anisotropy of about a factor 2. Quite recently a sizeable anisotropy 
of the thermal conductivity has been observed j?| which reaches 30% in the superconducting state. 



The anisotropy of the penetration depth is perhaps the most striking one since it is directly related to the nature of 
the superconducting state. Indeed || at zero temperature, the ratio between the superfluid densities A^~ 2 (0)/A~ 2 (0) 
is 2.3 , which is in good agreement with the anisotropy of the normal state resistivity and of the square of the plasma 
frequency g. Moreover we have also to take into account that in Y Ba,2Cu/±0% , which has twice as many chains as 
Y BaiCuzOi and an orthorhombic distorsion of 0.8% , the superfluid density anisotropy || is as high as 6.2 . This trend 
points toward the natural conclusion that the chains are conducting and that the superfluid density anisotropy is due 
to their direct contribution. Now we attribute the linear term at low temperature to nodes in the gap. Since in A^ 2 (T) 
this linear term is typically 3 times the linear term in A~ 2 (T), the nodes contribution is also strongly anisotropic. 
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This leads to the suspicion that chains have also something to do with the nodes. This is not so easy to include in 
the spin fluctuation picture because it deals only with the Cu02 planes. This gives rather a strong hint that we have 
to include chains explicitely in our modeling of the superconducting state in order to come up naturally with such 
an anisotropy. Another strong indication that chains play an important role in the superconducting properties comes 
from very recent penetration depth measurements (| which are partially in contradiction with earlier results [|l0| , and 
show a kink in the temperature dependence. This seems to indicate that a two component model is appropriate for 
the description of superconductivity and leads again to a two band model. Naturally the present overall experimental 
picture is still somewhat uncertain and future experiments will help to clarify the situation. 

These difficulties, as well as others, with the spin fluctuation mechanism have led us to propose recently a model 
[ pT[ which allows to account for all the above experimental facts. In this model pairing in the planes is due to an 
attractive interaction, while Coulomb repulsion induces in the chains an order parameter with opposite sign. It is well 
known that these plane and chain bands hybridize because electrons are physically allowed to jump from planes to 
chains. Due to the anticrossing produced by this hybridization, one obtains automatically an order parameter which 
changes sign on a single sheet of the Fermi surface and has nodes in the gap in the region of the Brillouin zone where 
the anticrossing takes place. In this sense our order parameter is d-wave like. However in contrast with the d x 2_ y 2 
order parameter, there is no symmetry breaking and we have a total of 8 nodes ( 4 on each sheet of the Fermi surface 
) instead of 4. Since the chains play an essential role in our model, it is meaningless to consider an approximate 
tetragonal symmetry. Rather we have to consider only the orthorhombic symmetry, under which our order parameter 
is completely invariant (no symmetry breaking), which is also the case for the d x 2_ y 2 order parameter . We can also 
say that on each sheet of the Fermi surface our order parameter looks qualitatively like a d+s order parameter. 

Since the order parameter of our model has nodes, it accounts naturally for all the experiments providing evidence 
for its change of sign and for the existence of low energy excitations. At the same time it provides also a simple 
explanation for the important anisotropy in the ab plane mentionned above, because of the important role of the 
chains. We have found [Q that the existence of two weakly coupled bands (plane and chain) in our model, as well 
as weak plane-chain scattering, leads to the possibility of a weak sensitivity of the critical temperature to impurities, 
in agreement with experiment. We find |]l3| that our model accounts quite well for the anisotropy of the penetration 
depth and for the absolute values. We are also able to reproduce fairly well the whole temperature dependence for 
both the a and the b directions found in Ref. juj, including the linear dependence at low temperature. 

Although our model is quite satisfactory for YBCO, in part because there are strong experimental indications that 
chains play a significant role, it appears for the same reason quite specific of YBCO and does not seem to be applicable 
to other compounds. However the natural tendency is to find a universal explanation for high T c superconductivity 
and hence a model applicable to all cuprate superconductors, since they all share fairly high T c and CuOi planes. On 
this basis our model sounds quite unsatisfactory. Let us merely remark that fortunately this universality hypothesis 
can be checked and hopefully will be checked experimentally. We note also that the available data do not seem to 
confirm completely this hypothesis, since N d,2~ x G e x CuO /± (NCCO) on one hand and YBCO (for the strong anisotropy 
discussed above) on the other do not fit in the picture. The case of NCCO will be discussed below. It would clearly not 
be satisfactory to find models specific to each of the different high T c compounds, but it is possible that, while sharing 
common features, they also differ in some respect. A related unappealing feature of our model is its complexity, since 
it is basically a two band model, not as nice as a one band model. Ideally one would certainly rather like to have a 
simple rather than a complicated model. However we already know that these compounds are chemically complicated. 
The same might unfortunately be true for their physics. 

Whatever the feeling about the above issues, there remains a more practical question. Other high T c compounds 
have been investigated experimentally and have shown features similar to the ones found in YBCO ( although the 
experimental picture is not as complete as in YBCO). Is it possible to extend the ideas of our model and find 
satisfactory explanations for them ? And more generally, following these ideas, what do we expect to find in other 
compounds ? It is the purpose of this paper to consider this question. As we mentionned our model for YBCO is 
basically a multiband model. Therefore we will in particular consider whether a multiband model is applicable to a 
given compound and if this can lead to reasonable explanations for experimental facts. We note that a multiband 
model applies trivially in all the compounds having more than one CuOi plane in the elementary cell. The usual 
trend is to consider that this leads to almost degenerate bands and to treat them as a single band. It is not completely 
obvious that this is correct, but anyway this will not be our main emphasis. We will rather be interested in looking if 
there are bands, analogous to the chain band of YBCO, that is not linked to the CuO^ planes, which play physically 
an important role. Naturally we do not intend to review the whole situation but only to make our best to discuss 
what we believe to be relevant points. Our conclusion is that we can extend our model in a natural way and we do 
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not find major incompatibility with experiments. Naturally this is a somewhat imprecise statement and we will be 
more specific below. We consider successively the cases of BSCCO, HBCCO, TBCCO, and finally LSCO and NCCO. 



II. BSCCO 



We turn first to the extension of our model to Bi2Sr2CaCu20%+s (BSCCO), which is probaby the second most 
investigated high T c superconductor after YBCO. This is not a very good compound for fundamental investigations 
since it is generally not stoechiometric but rather has often an oxygen "doping" 5 giving the maximum T c . This 
comes in addition to the incommensurate distorsion arising in the BiO planes. This is in contrast with YBCO where 
the nearly stoechiometric compound has been mostly investigated. We first note the high value of 2A/T*e (which is of 
order 7 or even higher, as obtained from tunnelling experiment |14| as well as from Raman scattering |L5|) gives an 



indication that a simple d-wave model does not properly describe BSCCO. This value is incompatible with a weak 
coupling d-wave model, and it is doubtful that strong coupling spin fluctuations calculations can lead to such a 
high value, because the spectrum for spin fluctuations goes to high energy and the coupling constant is of order unity. 
A related difficulty is the value of the slope of the penetration depth A(T) at low temperature p7| , |Tq] . In the n model 
of Xu et al. |l9| , which is an extension of the simple d-wave model allowing for a variable slope for the opening of the 
gap, A(0) = A fi for < 9 < fi^ 1 and A(0) = A for < 6 < tt/4 . The low temperature A(T) is given by [^9): 
[ A(T) - A(0) ] / A(0) = ( 2 ln2 / /j, A ) T . Taking the experimental value A(T) - A(0) « 10 A / K and T c w 90 K , 
this leads to /i Ao / T c = 1.5 10 -3 A(0) , where A(0) is measured in A . There is some uncertainty in A(0). If we take 
the assumed values of 2100 A Q and 2600 A@,we have either A / T c w 3.2 or [i A / T c w 4 . Taking the ( 
pretty unphysical ) lower limit fi = 4 /tt leads to 2Ao / T c = 5 or 6.2 , somewhat below experiment. The standard 
d-wave corresponds essentially to n = 2 and gives 2Ao / T c = 3.2 or 4 , barely compatible with weak coupling and in 
clear disagreement with experiment. It might naturally be that further experiment give a higher value for A(0) and 
a lower one for 2A/T C which would solve this difficulty. 

The natural way to extend our two-band model to BSCCO is to assume that the BiO planes lead to a metallic band, 
which would play the role of the chain-band in YBCO. The existence of this band at the Fermi level is controversial. 
There is some experimental evidence from photoemission that BiO planes are insulating because the BiO band is 
not seen [^0| , but these experiments are somewhat contradictory |^l[ . Similarly STM seems to indicate that the BiO 
plane at the surface is nonmetallic |l4|] . However both kind of experiments test only the surface of the sample on 
a few tens of A, that is basically one cell depth. Just because the surface necessarily implies some kind of lattice 
distorsion and/or consequently some modification in the electronic spectrum, the BiO planes may well be insulating 
at the surface while they are conducting in the bulk. It may also very well be that there is some oxygen depletion at 
the surface. An example of this kind of problem is found in YBCO ]22| . In this compound the bands corresponding 
to CuOi planes are seen in photoemission, while the CuO chain-derived band is not seen. This could lead again to 
the conclusion that this chain band is nonmetallic. However this chain band is seen by ACAR (which does not see the 
plane bands) which shows that this band is metallic. As indicated above this is confirmed completely independently 
by the in-plane anisotropy of the conductivity and mostly by the strong anisotropy of the superfluid density. Other 
problems arising with the interpretation of photoemission experiments are discussed in jy . 



The main reason to believe that the BiO planes are metallic comes from band structure calculations p3|]24| ] . Indeed 
these papers find that, in addition to the two nearly degenerate CuO^ plane-derived bands, there is at the Fermi level 
a BiO-derived band. The plane bands are rather similar to the ones found in YBCO, and the corresponding Fermi 
surface looks roughly as a circle centered at (± it / a , ± 7r / a ) if one uses the simplified Brillouin zone corresponding 
to tetragonal C11O2 planes. The BiO band gives a smaller Fermi surface centered at ( ± ir / a , ) and ( , ± ir / 
a ). Naturally one may wonder about the validity of band structure calculations. Clearly there are limitations in the 
precision of these complicated calculations. They are also sensitive to the precise location of the atoms. And LDA, 
which is commonly used to take into account correlations, is not an exact method for dealing with highly correlated 
systems. As a result mass renormalization due to interactions is not included and the calculations provide essentially 
the bare excitation spectrum (note however that mass renormalization due to Coulomb interaction seems to be rather 
small, from all the calculations dealing with correlation effects) . On the other hand the results from different authors 
are in fair agreement, which indicates at least that the results are reasonably well converged. Moreover in the case of 
YBCO, the results are in reasonable agreement with all experimental informations from ARPES, ACAR and dHvA 
p2fl , including a small electron pocket centered at (± n / a , ± n / a) , which is apparently not of physical importance. 
It is therefore difficult to ignore these results and one has to consider seriously the possibility that their conclusion is 
correct. 



3 



In the case of BSCCO the BiO band is mostly above the Fermi level, but the bottom of the band is located around 
0.5eV below the Fermi level. It is likely that the precision of the method is better than this, and the BiO band will 
not be pushed above the Fermi level by some imprecision in the calculations. Therefore the existence of this BiO 
Fermi surface seems rather secure. On the other hand the calculations show a second BiO band which comes quite 
close to the Fermi level from above (within O.leV or less), but does not cross it. It is not obvious that the precision is 
good enough to make sure that there is not a second BiO sheet of the Fermi surface. Experimentally, as mentionncd 
above, ARPES does not see the BiO band, but it sees the CuOn bands ( or at least one of them, but the failure to 
see two bands might be due to their splitting being very small ) . The Fermi surface corresponding to this band is in 
good agreement with band structure calculations, which gives again confidence in these calculations. On the other 
hand ACAR [£5| gives results in fair agreement with band structure calculations and sees the BiO piece of the Fermi 
surface, although its size is somewhat larger than the theoretical prediction and the separation between BiO and 
CUO2 Fermi surfaces larger than predicted. Recent experiments pq] on the mid-infrared reflectivity in Bi 2 Sr 2 CuOe 
have also been interpreted in terms of transitions between BiO and CuO-i type bands. As mentionned above the 
discrepancy with ARPES might be due to the fact that the surface structure is different from the bulk. A similar 
problem may be seen in SY2-RUO4, a layered perovskite becoming superconducting around IK. Here there is some 

] and dHvA (28) experiments, and it has been suggested (2^]that the origin might 
be surface relaxation of the oxygen position. 

Actually calling this electron pocket the BiO piece of the Fermi surface is a misnomer. Indeed the corresponding 
electronic states are J2^]24|] strong hybrid of states describing electrons moving respectively in BiO planes and in CuO-i 
planes. The situation is quite similar to the one found in YBCO where we start with electrons belonging either to 
planes or chains, and which are hybridized due to the hopping of an electron from plane to chain, through the apical 
oxygen. In the same way we can start in BSCCO with "non-interacting" bands, as considered explicitely by Massida 
et al. pi) , corresponding to electrons moving either in BiO planes or in CuOi planes. Then hopping through the 
apical oxygen produces an hybridization of these bands Jm). The net result is that (24) the BiO planes contribution 
to the total density of states at the Fermi level is almost half of the one of the C11O2 planes, which shows clearly the 
importance of BiO planes. 

It is possible to produce semiquantitatively the results of band structure calculations with a simple tight-binding 
model. Let us indeed consider that BSCCO is build by stacking up sets made of two C11O2 planes and two BiO planes. 
Motion in the C11O2 planes and in the BiO planes is described by the dispersion relations €k and e' k respectively. 
Moreover we allow hopping toward the nearest plane. We call t the hopping term between a Cu0 2 plane and the 
nearest BiO plane, t' the hopping term between the two BiO planes and t" the one between the two C11O2 planes, 
assuming for simplicity that they are wavevector independent. The corresponding Hamiltonian is: 

#0 = E £fcC+i,n C fc,l,« + J2 £k4 2 . n Cka,n + E *" (Cfe,l,n C M,« + h -C.) 
n n n 

+ E £'kdt thn d k ,l,n + E £ ' 'fcdfc,2,rA2,« + E *' ( rf fc,l,„+l4,2,n + h.C.) (1) 
n n n 

+ E t ( C k.2,n d k,l,n + h.C.) + E t [c£ fl>n dk,2,n + h.C.) 
n n 

where c k and dt are creation operators in the C11O2 planes bands and in the BiO planes bands respectively. The 
indices 1 and 2 number the Cu0 2 planes and BiO planes, while the index n numbers the stacks. Actually we can 
simplify further by taking into account that BSCCO is very anisotropic which implies that some hopping terms are 
quite small. From band structure calculations Q it is clear that hopping between C11O2 planes and between BiO 
planes are small, while hopping between neighboring Cu0 2 plane and BiO plane are somewhat stronger, leading to 
rather isolated Cu0 2 - BiO stacks. For simplicity we assume for example hopping between BiO planes negligible, 
that is t' = 0, and t" ^ . But we could take the opposite assumption t' ^ 0, and t" = and reach the same 
qualitative conclusions. The important point is that the sets made of two Cu02 planes and two BiO planes are now 
independent and we have a two-dimensional problem. The Hamiltonian can then be separated in two independent 
parts H± corresponding to even and odd combinations. Introducing the operators c± and d± by ci,2 = ( c+ ± c_ 
)/y/2 and di, 2 = ( d+ ± d- )/y/2 : 

H± = ^(£ fe ± f )cl ± c K± + £'kdt ± d k ,± + Y,t{c+ ± d k ,± + h.c.) (2) 

Each of these Hamiltonian is analogous to the two-dimensional model |ll|] we have considered for YBCO. For each of 
them the Cu0 2 band and the BiO band hybridize. Since band structure calculations show that one of the BiO derived 



disagreement between ARPES |27 
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band is pushed above the Fermi level by this coupling, we consider only the other Hamiltonian, where hybridization 
leads to an anticrossing of the Fermi lines. In order to obtain a Fermi surface in reasonable agreement with band 
structure calculations, we take for the unhybridized CuOi band the rather standard dispersion relation : - t" = 
- 2£o ( cos^a) + cos(fcj,a) ) + 2to cos(fc x a) cos(fcj,a) - /i with to — 0.33 eV, fj, = - 0.46 eV (the axes are along the 
CuO bonds). On the other hand the unhybridized BiO band can be described by a nearest neighbour tight binding 
approximation (with BiO bonds rotated by 45° with respect to the CuO bonds), which leads to e' k = 4t' cos(k x a) 
cos(k y a) - [i! , with t' = 0.3 eV and fjf = - 0.6 eV in order to reproduce band structure calculations. The result for 



the Fermi lines is shown in Fig.l for t = 0.1 eV . This figure is very similar to the Fermi surfaces found in |!3 24|. We 
can then take over our model @ for YBCO and assume that there is a repulsive pairing interaction between C11O2 
plane and BiO plane, leading to an order parameter with opposite sign in the CuOi plane and BiO plane. Then 
because of hybridization we have nodes in the anticrossing regions. To give an example we assume the regime of well 
separated bands, where the hopping term t is large compared to the gaps (actually it is not clear that this condition 
is satisfied because the maximum gap is quite large in BSCCO as we have seen above, but this should not lead to 
qualitative changes). In this case, assuming for simplicity isotropic order parameters A in the CuOi plane and A' in 
the BiO plane, the order parameter is given jl2| by A& = ( A e' k + A' ek ) /( £fc + e' k ) when hybridization is taken 
into account. In particular the locations of the nodes is not so sensitive to the ratio A/A' because the Fermi lines for 
unhybridized C11O2 plane and BiO plane cross essentially at right angle. We show in Fig.l the locations of the nodes 
for A/A' = - 1 . 

We can then try to calculate the temperature dependence of the penetration depth, just as we have done for YBCO 
p3| . However the results depend sensitively on a number of ingredients |l3| ] and the situation is here much more 
uncertain than in YBCO, in particular with respect to the band structure and the parameters describing the BiO 
derived bands. Therefore we will only consider the low temperature behaviour where the number of hypotheses can be 
somewhat more limited. We proceed as in |]l3| and treat the hopping term t as a small quantity so that hybridization 
modifies the dispersion relations only in the vicinity of the crossing of the unhybridized Fermi lines. This leads to : 

A-(0) - X-\T) = —^-—^Av (3) 
7r U 2 c(A|A'|) 1/2 

where A and A' are the order parameters in the unhybridized CuOi and BiO planes respectively (which we have 
assumed to be k-independent for simplicity) and c = 14.6 A is the thickness of the stack made of two C11O2 planes 
and two BiO planes. The dimensionless constant Aq is related to the quasiparticle velocities at the crossing point 
of the unhybridized Fermi lines. We assume for simplicity that these velocities are orthogonal (this is essentially the 
case in Fig.l). In this case Aq is given by Aq = (A 2 v g + A' 2 v c )/(A |A'| vb i>c) where vc and vb are the velocity 
of the unhybridized CuOi and BiO bands respectively. This quantity has a minimum equal to 2 when A vb = |A'| 
Vc and we take for simplicity this worst case for our evaluation. Expressing numerically all known physical constants, 
this can be rewritten in the low T regime as : 

X(T)-X(0) = 2.8m^^-^t (4) 

where the all penetration lengths are now expressed in units of 1000 A and the hopping parameter t is in eV. We want 
now to compare this result with experiments on BSCCO which gives a slope of 10 A/K. We take for example 

A/ A' = - 1 (actually we would expect A' to be somewhat less than A, which would be make our case easier), and 
A(0) = 2000 A from experiment (again a rather worst case hypothesis). If we use A/T c w 3 as discussed above, we 
obtain agreement with experiment for t = 0.12 eV . This result is quite similar to what we have found for YBCO. It 
is small enough to make our hybridization approach consistent. This is quite satisfactory since we have basically no 
adjustable parameters. On the other hand it is not in contradiction with the strong anisotropy and the large resistivity 
well known in BSCCO, since this can be attributed entirely to the very small hopping between either the BiO planes 
or the CuOi planes or both, as we have already discussed. It may seem rather strange that we can account for the 
rather strong slope found experimentally by a small hybridization. Roughly speaking this is due to the fact that we 
have a total of 16 nodes in the Brillouin zone, whereas the standard d-wave model has only 4 nodes. 

Our model differs clearly from standard d-wave, not only by the number of nodes but also by their locations. This 
is in contrast with YBCO where the location of the crossing between plane and chain bands is not so far from the 
45° direction of the standard d-wave model. Here, in Fig.l , the directions of the nodes are respectively 18° and 
22° (and all the other ones obtained by tetragonal symmetry), with the origin taken at (tt, 7t) as it is usually done. 
This is quite different from the 45° of the d-wave model. Naturally Fig.l is just an example, and it is likely that 
the situation in actual BSCCO is somewhat different because band structure ( as suggested by ACAR experiments) 
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and order parameter are not the ones we have chosen. However we do not expect the location of the nodes to 
change qualitatively. Therefore experiments designed to ascertain the position of the nodes are quite important and 
we consider now briefly the evidence for this position in BSCCO (we will also discuss this point below for HBCO). 
At first photoemission experiments were interpreted |3(| as showing nodes at 35°. Taking together into account 
experimental and theoretical uncertainties, we would not consider this as being in clear contradiction with our model. 
However these experiments have then been reinterpreted [fjlf as giving the nodes at 45° which agrees with d-wave and 
disagrees with our model. Anyway, as we discussed above, it is not clear at all that photoemission provides a faithfull 
picture of what is happening in the bulk. And it sees only a Cu0 2 derived band, and no BiO band. So the fact that 
our model disagrees with photoemission is not so surprising. 



On the other hand in-plane tunneling experiments give an indication on the location of the nodes. They have 
shown a clear gap anisotropy 32 1 with a minimum gap along the CuO bond direction in contrast with photoemission 
experiments. Naturally, just as photoemission, this technique is limited in angular resolution and can not by itself 
display nodes, but the obvious conclusion is that the nodes of the gap are not so far from the CuO bond direction, 
if there are nodes at all in BSCCO. This is not in agreement with d-wave and more in favor of our model. Finally 
half-integer flux quantum effect has been observed in tricrystals of BSCCO j33|, providing rather convincing evidence 
for a change of sign of the order parameter on the Fermi surface. However this experiment does not provide a unique 
answer for the location of the nodes. For example, although the geometry is essentially the same as the one used for 
YBCO and is therefore consistent with d-wave, a g-wave order parameter A(#) « cos(46*) is also a possible solution. 
More generally, since the Josephson junctions ( and in particular the scattering at the interfaces ) are not controlled 
in the tricrystals, it seems quite difficult to be secure about the position of the nodes from this kind of experiments. 
In conclusion it is fair to say that the problem of locating the nodes in BSCCO is not at all settled. At the same time 
it would clearly be a crucial piece of information to know in a reliable way this location since this would definitely 
allow to narrow the number of possible models and be another step toward identifying the mechanism of high T c 
superconductivity. 



III. HBCCO 



Let us consider now the case of the Hg compounds. Band structure calculations give a more complex situation than 
for BSCCO. Indeed the compound HgBa 2 Cu04 + s (Hg-1201) is found to have only a Cu0 2 derived band crossing 
the Fermi level [^4]-p7|, although there is a HgO derived band which is coming quite near ( typically 0.1 eV) to cross 
also the Fermi level from above. It is not clear that the precision of the calculations is enough to be certain about 
it. Moreover the role of oxygen doping, which is held responsible for the metallic and superconducting properties of 
this compound [B4[, is also difficult to take into account and it seems to be quite important |p8|. Therefore it will 



be quite interesting to see what experiment will tell in this respect. The present tunneling data 39 are more in 
favor of a standard BCS density of states. This is coherent with our model since, with a single band, we expect only 
anisotropic s-wave pairing. However a two-band model is still an open possibility in this compound. An indication 
in this direction seems to be that, for Hg nuclei, T\ has essentially the same value at T c for Hg-1201 and Hg-1212, 
showing a similarity between these two compounds. 



Indeed calculations in HgBa 2 CaCu 2 6+s (Hg-1212) and HgBa 2 Ca 2 Cu 3 8+s (Hg-1223) find |§|3|@ the HgO 
derived band crossing also the Fermi level, in addition to the Cu0 2 bands, which makes these compounds "self-doped" . 
The situation is actually quite similar to the one found in BSCCO. The Fermi surface displays a small electron pocket 
centered at the X point. Because of hybridization between HgO and Cu0 2 bands, the electronic states on this sheet 
of the Fermi surface have a mixed HgO - Cu0 2 character, as well as those of the nearby Cu0 2 bands in the vicinity 
of the X point. An anticrossing is clearly seen. Therefore the necessary ingredients for our model are present, just as 
in BSCCO, and they can lead in the same way to nodes in the gap. Experimentally Raman scattering on Hg-1212 
jyjgives evidence for nodes in the gap, but it shows also that some nodes (if not all) are away from the 45° direction. 
This is in disagreement with d-wave and in agreement with our model. The low temperature dependence of the 
in-plane penetration depth in Hg-1223 (42jseems to be linear, which would indicate the existence of nodes in the 
gap. The available tunneling |^3| data are also suggestive of nodes. To summarize the present situation in these Hg 
compounds is coherent with our model. 



G 



IV. TBCCO 



We consider finally the case of the Tl compounds. Here again band structure calculations give results quite 
similar to the ones found for Hg compounds. There is a TIO band which crosses |44|-|47j| the Fermi level for 
TliBa2Ca n -\Cu n O/±j r 2n with n=l,2,3 and 4, giving rise to a small electron pocket for the Fermi surface. How- 
ever for TliBa2CuO§ there is some disagreement between various calculations Jf4],(l7| and this result seems to be 
within the error bars. There are differences with Hg compounds. First the electron pocket is found around the T 
point. The states corresponding to this part of the Fermi surface are fairly delocalized throughout the cell with 
contributions not only from Tl and O atoms, but also from the apical O and even the d{z 2 ) orbital of Cu atoms ( as a 
result, for TliBa2CuO§ E7] |, this piece of the Fermi surface has a more pronounced 3-dimensional character ). Then 
there is also the possible existence of hole pockets |^J45| , derived from flat CuOi bands. Finally since one has always 
the problem of dealing with oxygen doping, compounded with the experimental uncertainty on the atomic positions, 
the band structure situation is quite complicated. 

We can again apply our model as in the case of BSCCO and say that Coulomb repulsion will lead to an order 
parameter with opposite signs on the sheets of the Fermi surface linked to the CuOi band and to the TIO band. 
This might explain the observation |48[ ] of spontaneously generated half quantum flux in TI2BCL2C11O6+S tricrystals. 
Nevertheless a major difference in the results of band structure calculations as compared with the Bi and Hg compounds 
is the lack of crossing of the C11O2 and the TIO bands. Therefore in our model we would not expect nodes to appear. 
However it might very well be that the experimental situation is in this respect slightly different from the one 
predicted by band structure calculations |44|j47|] and that there is actually some crossing. This possibility has also 
been considered in the context of the inter-layer pairing mechanism ]49| ] . As we have mentionned there is enough 
uncertainty in the band structure calculations to allow for this possibility. Clearly direct information from ARPES, 
ACAR,dHvA experiments would be quite useful in order to settle this question. In our model these crossings would 
lead to nodes in the gap, in a way analogous to what we have found in BSCCO and HBCO. Actually available Raman 
scattering experiments favor the existence of nodes in a way very analogous to HBCO since they give, in B± g 
symmetry, a response linear in frequency at low u>, in the superconducting phase. 

V. NMR 

Since our model requires two kind of metallic bands in the Bi, Hg and Tl compounds, it appears important to find 
experimentally if we are indeed in this situation or not. This is actually rather difficult. We have already discussed this 
problem for BSCCO and found that there is no good experimental evidence against metallic BiO planes. Naturally 
one would rather like to find experiments proving definitely that they are indeed metallic, but we are not aware of 
any. Surprisingly the (quasi) quadratic structure of these compounds is in this respect a disadvantage, since in the 
case of YBCO the most reliable evidences of the metallic character of the chains come from the anisotropic properties 
linked to the orthorhombic structure. The difficulty comes from the microscopic nature of the information we are 
looking for, since for any macroscopic property we will not be able to know which one of the two bands is responsible 
for it. 

Nevertheless it would seem that the NMR is an effective experiment in order to obtain this information since it is a 
local probe. The case of YBCO is a good example of this pSl)] , although it seems to be somewhat forgotten. Indeed the 
superconducting transition is clearly seen as a strong drop both in the Knight shift K and in the relaxation rate I /T± 
of Cu(l) nuclei of the CuO chains, and the behaviour below T c is qualitatively similar to what is observed for Cu(2) 
in the CUO2 planes. This shows that there is on this Cu(I) site a non vanishing amplitude for the wavefunction of the 
quasiparticles subject to the superconducting condensation. Moreover the quantitative temperature dependence for 
Cu(I) below T c is clearly different from the one observed for the planes Cu(2), both for K and 1/T% . This shows that 
what is seen by Cu(l) nuclei is not merely a transferred part of the wavefunction of quasiparticles located dominantly 
in the C11O2 planes, since in this case we would expect the temperature dependence to be the same for Cu(I) and 
Cu(2). Therefore the NMR favors a two-band picture for YBCO. This picture is also supported by the fact that the 
magnitude of at T c are quite similar for Cu(I) and Cu(2), and the same is true for the change in Knight shift 
between T = and T c . 

Let us consider the situation for Hg compounds. 199 Hg NMR experiments have been performed recently [Qon 
HgBa,2CuOi+s . Below T c they display a strong drop both in the Knight shift K and in the relaxation rate \JT\ , 
so the Hg nuclei see the superconducting electrons. Moreover for 199 Big , T\T is essentially constant in the normal 
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state ( the Knight shift is also constant with a "Korringa ratio" rather near what is expected from a Fermi liquid ), 
whereas (TiT) -1 increases markedly with decreasing temperature for 63 Cu. This is in favor of a two-band picture, 
just as in YBCO. One may argue that antiferromagnetic fluctuations are responsible for the increase in ( 63 TiT) _1 
and that these fluctuations completely cancel for symmetry reason at the Hg site, so that Hg has only a Korringa 
behaviour. However this implies an antiferromagnetic coupling between different CuOi layers which is surprisingly 
strong ( taking in particular into account the disorder induced by oxygen doping in the Hg planes ). A puzzling 
ingredient to this question is provided by the NMR on HgBa2CaCu20e+s ]53|] where the temperature dependence 
of (TiT) -1 for 63 Cu and 199 Hg is nearly the same. Clearly further experiments are needed to elucidate this point. 
Finally one has 199 Ti w 10 2 . 63 Ti , but it is much more difficult to obtain firm conclusions from a comparison of the 
magnitude of 63 Ti and 199 Ti since we deal with different nuclei, in contrast to YBCO. 



We turn then to the Tl compounds. Again in T^BaiCuO^+s |54 55|thc Knight shift as well as the relaxation rate 
of 205 Tl display below T c a strong drop, which shows that the wavefunction of superconducting electrons extends in a 
quite sizeable way to Tl nuclei. One may wonder, as for the Hg compounds, if this is not due to a transferred hypcrfine 
interaction of the Hg nuclei with electrons on the Cu site of the nearby C11O2 planes through the apical oxygen. This 
possibility remains to be investigated in details. However, just as in the Mila and Rice phenomenological analysis 
p6| for YBCO, one expects this hyperfine coupling to occur through the contact interaction due to Tl-s orbitals and 
therefore to be isotropic. The other contributions should be much smaller. However the spin part of the Knight 
shift has a sizeable anisotropy ( 30% in a T c = 110K sample |55| ), which is in contradiction with this hypothesis. 
Moreover, although it is difficult to obtain firm conclusions from magnitude considerations, it is rather striking to 
see that the 63 Cu and 2Q5 Tl relaxation rates are nearly the same ( 63 (l/TiT) w 20 sec^R- 1 and 205 (l/2iT) « 8 
sec -1 ir _1 in the normal state |5^j5^] ). Such a short 205 Ti is difficult to reconcile with a transferred hyperfine field. 
The natural interpretation of these TZ NMR results is rather a confirmation of the two band picture given by band 
structure calculations. 



Although more work is needed in order to reach a definite conclusion, we believe that the present NMR experimental 
evidence is leaning toward a two-band picture. However, even if we consider the single band hypothesis (assuming 
that only the CuO^ band is relevant), the NMR shows unambiguously that the electronic wavefunction extends in 
a significant way to the Hg or Tl sites. Just as in our two band model, this provides a possible physical origin for 
an important Coulomb repulsion within this single band. This repulsion is likely to be quite anisotropic, since the 
probability of finding an electron on the Hg or Tl site will clearly change when one moves over the Fermi surface. 
This anisotropic repulsion might be responsible for a change of sign of the order parameter within this single band. 
Therefore in any case the NMR provides a clear hint that the Hg or Tl sites can not be safely omitted in a proper 
description of the electronic properties of the Hg or Tl compounds. 



VI. NCCO AND LSCO 



Let us finally consider the high Tc compounds, such as La,2- x Sr x CuO '4 or Nd2— x Ce x CuO& , which display only 
CuOi planes without any other possible metallic component and hence provide a check for our physical picture. 
Indeed in this case we are naturally led to a single band and naturally band structure calculations lead in these 
compounds to a single C11O2 derived band |57|. In the framework of our model we do not expect any nodes in the 
gap. This leads for these compounds to s-wave pairing, even if it is quite possible that the order parameter is fairly 
anisotropic. The experimental evidence in Nd,2- x Ce x CuOi is presently clearly in favor of s-wave pairing. Indeed the 
low temperature dependence of the penetration depth follows the standard weak coupling behaviour 58 1 . Tunneling 
gives also a fairly standard BCS shape, and moreover the Eliashberg function extracted from these tunneling 
data leads to a T c in good agreement with experiment, providing evidence for phonon mediated attractive interaction. 
Finally recent Raman scattering experiments [ |60| go also in the direction of s-wave pairing since the results can be 
explained by an almost isotropic gap. Therefore NCCO seems in contradiction with the spin fluctuation mechanism. 
The fact that it is an "electron doped" compound, in contrast with all the other "hole doped" compounds is a rather 
unconvincing explanation for this discrepancy. Indeed one starts with a Fermi liquid with strong spin fluctuations 
and it is not clear at all why a small change in doping, allowing to go from the hole doped side to the electron doped 
side, would kill these spin fluctuations or reduce strongly their coupling to the electrons, while at the same time a 
sizeable coupling to phonons would appear. 



On the other hand the experimental situation is unclear in La,2- x Sr x CuOi with some experiments in favor of s-wave 
and others in favor of nodes. To our knowledge there is no phase sensitive experiments, nor precise measurements 
of the low temperature behaviour of the penetration depth. Recent Raman scattering experiments have found 
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results consistent with d-wave symmetry. However these results are not precise enough to prove convincingly the 
existence of nodes in the gap and some of their features do not agree with d-wave, but they display a clear anisotropy 
in the gap. Anyway, whatever the model, an experimental conclusion giving nodes in LSCO and no nodes in NCCO 
would be quite puzzling and we would rather hope for our understanding of high T c superconductivity that both 
LSCO and NCCO fall in the same class, whatever it may be. 



VII. CONCLUSION 



In this paper we have proposed an extension to other high T c compounds of a model introduced earlier for YBCO. 
This extension relies on the fact that a number of these compounds are believed to be " self-doped" , with a part of 
the compound giving holes to CuC>2 planes. Explicitely this doping part corresponds to the BiO, HgO, TIO planes in 
BSCCO, HBCCO, TBCCO respectively, just as it is the CuO chains in YBCO. This doping part is then considered 
as inert or unimportant in many models. We propose instead that it is metallic, which leads naturally to a multiband 
model. This picture is supported by band structure calculations. Naturally these calculations can not make sure 
that these doping bands are indeed metallic, but they make this possibility quite plausible. The final answer on this 
point can only be obtained from experiment. We can then take over the ingredients of our model in YBCO, and 
assume a repulsive pairing interaction between these doping bands and the CuO?. bands. This leads to an opposite 
sign for the order parameter on these bands and to nodes whenever the Fermi surfaces of these bands cross. We have 
considered in particular the case of BSCCO and shown that this model accounts with reasonable parameters for the 
low temperature dependence of the penetration depth. An interesting feature of BSCCO is that the nodes arc not 
located near the 45° direction found in the d-wave model. This is in contrast with YBCO where the location is not 
far from 45°. Therefore the reliable experimental determination of the node locations is an important test for our 
model (we have discussed why ARPES does not provide a final answer in our view). We have then discussed the 
cases of HBCCO and TBCCO which have not been explored as much as BSCCO up to now. The situation in these 
compounds is somewhat similar to BSCCO with respect to our model. We have also considered the experimental 
evidence given by NMR and found that there are good indications for a metallic character of the doping bands in 
these last two compounds, although more experimental work is needed in order to get a definite anwer. Finally our 
model does not give nodes for NCCO and LSCO. This is in agreement with the present experimental situation for the 
first compound, but the evidence is quite ambiguous for the last one. We hope that experimental progress will allow 
fairly soon to confirm or invalidate our description. However whatever the answer we definitely need a consistent 
picture of all high T c superconductors, in agreement with all experimental data. 
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FIG. 1. Fermi surface in the Brillouin zone. The filled circles give the positions of the nodes of the gap. 
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